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ABSTRACT 

Dynamical mass measurements to date have allowed determinations of the mass M and the distance 
D of a number of nearby supermassive black holes. In the case of Sgr A*, these measurements are 
limited by a strong correlation between the mass and distance scaling roughly as M ~ D 2 . Future very- 
long baseline interferometric (VLBI) observations will image a bright and narrow ring surrounding the 
shadow of a supermassive black hole, if its accretion flow is optically thin. In this paper, we explore the 
prospects of reducing the correlation between mass and distance with the combination of dynamical 
measurements and VLBI imaging of the ring of Sgr A*. We estimate the signal to noise ratio of near- 
future VLBI arrays that consist of five to six stations, and we simulate measurements of the mass 
and distance of Sgr A* using the expected size of the ring image and existing stellar ephemerides. 
We demonstrate that, in this best-case scenario, VLBI observations at 1 mm can improve the error 
on the mass by a factor of about two compared to the results from the monitoring of stellar orbits 
alone. We identify the additional sources of uncertainty that such imaging observations have to take 
into account. In addition, we calculate the angular diameters of the bright rings of other nearby 
supermassive black holes and identify the optimal targets besides Sgr A* that could be imaged by a 
ground-based VLBI array or future space- VLBI missions allowing for refined mass measurements. 
Subject headings: accretion, accretion disks — black hole physics — galaxy: center — galaxies: nuclei 
— gravitational lensing: strong 



1. INTRODUCTION 

Sgr A*, the supermassive black hole at the center of 
our galaxy, has been observed for several decades. Mon- 
itoring stars orbiting around Sgr A* has led to mea- 
surements of its mass and distance (Ghez et al. 2008; 
Gillessen et al. 2009). However, these measurements of 
mass and distance are strongly correlated. For purely as- 
trometric measurements, mass and distance are related 
as M ~ D 3 , while for measurements of radial velocities 
mass and distance are related as M ~ D . For combined 
data sets, the correlation between mass and distance be- 
haves roughly as M ~ D 2 (Ghez et al. 2008; Gillessen et 
al. 2009). This correlation between mass and distance 
constitutes a major source of uncertainty in our under- 
standing of the properties of Sgr A* . Likewise, dynamical 
measurements of the masses of a number of nearby super- 
massive black holes have been obtained with often much 
greater uncertainties (see, e.g., Giiltekin et al. 2009). 

Another technique, VLBI, aims to image Sgr A* di- 
rectly. Recent VLBI observations with an array consist- 
ing of the Submillimeter Telescope Observatory (SMTO) 
in Arizona, the James Clerk Maxwell Telescope (JCMT) 
on Mauna Kea, and several of the dishes of the Com- 
bined Array for Research in Millimeter- wave Astronomy 
(CARMA) in California resolved Sgr A* on scales com- 
parable to its event horizon and identified sub-horizon 
size structures (Doeleman et al. 2008; Fish et al. 2011). 
Images of accretion flows around black holes have the 
shadow of the compact object imprinted on them, which 
depends uniquely on its mass, spin, and inclination (e.g., 
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Falcke et al. 2000) as well as on possible deviations from 
the Kerr metric (Johannsen & Psaltis 2010). Based on 
such images and assuming the mass and distance ob- 
tained from the monitoring of stellar orbits, these VLBI 
observations inferred constraints on the inclination and 
spin of Sgr A* (Broderick et al. 2009, 2011) and placed 
limits on potential non-Kerr signatures (Broderick et al. 
2012). 

In addition to the shadow, images of optically thin 
accretion flows around black holes carry a characteris- 
tic signature in the form of a bright ring (Johannsen 
& Psaltis 2010), which we refer to as the photon ring. 
Light rays that approach the event horizon closely orbit 
around the black hole many times before they are de- 
tected by a distant observer, resulting in a bright ring- 
due to their long optical path length through the ac- 
cretion flow. The flux of such photons can account for 
a significant fraction of the total disk flux and produce 
higher order images (Cunningham 1976; Laor, Netzer, & 
Piran 1990; Viergutz 1993; Bao, Hadrava, & Ostgaard 
1994; Cadez, Fanton, & Calvani 1998; Agol & Krolik 
2000; Bcckwith & Done 2005). These photon rings are 
clearly visible in all time-dependent general-relativistic 
simulations of accretion flows that have been reported to 
date (Moscibrodzka et al. 2009; Dexter, Agol, & Fragile 
2009; Shcherbakov & Penna 2010). 

Johannsen & Psaltis (2010) showed that a measure- 
ment of the ring diameter measures the ratio M/D for 
the black hole, independent of its spin or deviation from 
the Kerr metric. Therefore, combining such a measure- 
ment with the observations of stars around Sgr A* can 
reduce the correlation between mass and distance. 
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In this paper, we explore the ability of this approach 
to refine the mass and distance measurements of Sgr A*. 
We estimate the precision with which a thermal noise- 
limited VLBI array can infer the diameter of the ring of 
Sgr A* and use a Bayesian technique to simulate mea- 
surements of the mass and diameter of Sgr A* in conjunc- 
tion with parameters inferred from the existing data of 
the orbits of stars at comparable wavelengths. We show 
that, in this best-case scenario, the correlation between 
mass and distance is reduced significantly. In addition, 
we argue that the accretion flows of other nearby super- 
massive black holes are optically thin, allowing for VLBI 
observations of their respective photon rings. We assess 
the prospects of using this technique to infer the masses 
of these sources. 

2. MEASURING THE PHOTON RING OF SGR A* 

The properties of photon rings are practically indepen- 
dent of the specific flow geometry and remain constant 
even if the accretion flow itself is highly variable (Jo- 
hannsen & Psaltis 2010). The relative brightness as well 
as the constancy of these rings make them ideal targets 
for VLBI-imaging observations. 

For a Kerr black hole with mass M, the shape of a 
given photon ring has a diameter of 

L ~ 10.4 r g , (1) 

which remains practically constant for all values of the 
spin and disk inclination (Johannsen & Psaltis 2010). In 
this expression, 

r 9 = — (2) 

is the gravitational radius, and G and c are the gravita- 
tional constant and the speed of light, respectively. 

The angular diameter 6> r ing of the diameter of the pho- 
ton ring of a black hole is given by the ratio of its diam- 
eter and distance, 

#ring = yj- (3) 

Assuming the current mass and distance measurements 
of Sgr A*, M = 4.3 x 10 6 M Q and D = 8.3 kpc 
(Gillessen et al. 2009), the photon ring has an angular 
diameter of 

9 — 53 ^iarcsec. (4) 

Radio interferometers are limited by their intrinsic 
resolution as well as by interstellar scattering. In or- 
der to identify the range of wavelengths within which 
VLBI measurements of the photon ring of Sgr A* are 
resolution-limited, we compare the blurring effects of in- 
terstellar scattering with the resolution of an interferom- 
eter. In Figure Q] we plot the minimum size of resolvable 
structures on the image of Sgr A* using the interstellar 
scattering law of Bower et al. (2006). We also estimate 
(dashed line) the resolution of a radio interferometer at 
a given wavelength A by the expression 
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6> ros = kX/d 
with k = 1 and a diameter d = 10 4 km 



rth 



(5) 
which 



is comparable to the baseline length between the JCMT 
on Hawaii and the South Pole Telescope. This yields 
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Fig. 1. — Solid lines: Degree of the blurring of structures on the 
image of Sgr A*, inferred using the interstellar scattering law of 
Bower et al. (2006); the two lines correspond to the major and 
minor axes of the scattering ellipse. Long dashed line: Resolution 
estimate of a radio interferometer, taken as 8 Tes = kX/d with k = 1 
and d = 10 4 km. Dotted lines: Angular diameter of Sgr A* corre- 
sponding to length scales of 10r 9 and lr g , respectively. Imaging 
the photon ring of Sgr A* with VLBI is possible at wavelengths 
A < 1mm. 

/iarcsec and A mm is the observed wavelength in mm. Dot- 
ted lines mark the angular diameters corresponding to 
length scales of 10r g and lr g at the assumed distance 
of Sgr A*. As can be seen from this figure, at sub- 
mm wavelengths interstellar scattering becomes negli- 
gible and measurements are limited by the resolution. 
Therefore, a measurement of the diameter of the photon 
ring of Sgr A* will require VLBI observations at wave- 
lengths A < 1mm. 

In the following, we estimate the observed flux of the 
photon ring of Sgr A* assuming a Schwarzschild black 
hole and employing a model of a geometrically thin 
advection-dominated accretion flow (ADAF; Narayan 
& Yi 1994, 1995; Narayan, Yi, & Mahadevan 1995). 
ADAFs model the accreting gas around a supermassive 
black hole as a quasi-spherically symmetric plasma con- 
sisting of thermal electrons and ions at different tempera- 
tures (Narayan & Yi 1995a, b) as well as of non-thermal 
electrons (Mahadevan 1999; Ozel, Psaltis, & Narayan 
2000; Yuan, Quataert, & Narayan 2003). Such an accre- 
tion flow is allowed to cool through comptonization and 
the emission of bremsstrahlung and synchrotron radia- 
tion, with the latter generating the predominant contri- 
bution to the observed spectrum of Sgr A* at radio and 
sub-mm wavelengths (Narayan et al. 1995). 

For our estimate, we follow Broderick et al. (2009) 
and assume an ADAF model with a density of thermal 
electrons 



n e (r) = n e0 



electron temperature 
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where 9 r 



is the resolution in 



T e (r)=T eQ - 



(6) 
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and magnetic field 



B 2 n , 2GM m v 
8tt h 12r 1 



(8) 



where /3 = 10. For the coefficients we use (Broderick et 
al. 2011) 

3 x 10 7 cm" 3 (9) 



and 



T e0 = 1.7 x 10 11 K 



(10) 

These coefficients lead to predictions of the spectrum, 
polarization, and image size for Sgr A* that are in agree- 
ment with all current observations. 

We assume that all the emission at mm-wavelengths 
is due to thermal synchrotron radiation (see Narayan et 
al. 1995). While there may be small contributions of 
synchrotron emission from non-thermal electrons (e.g., 
Mahadevan 1998; Ozel et al. 2000; Yuan et al. 2003) or 
of other types of radiation at these wavelengths (such as 
jets; e.g., Falcke et al. 1993), this assumption affects our 
analysis only marginally. In practice, the measured total 
flux incorporates all such contributions. 

Following Dolence et al. (2009), we write the syn- 
chrotron emissivity as 

3» * ff^ia^ 1 ' 2 + 2 11/12 X 1/6 ) 2 exp(-AV3), 
3c-K 2 (l/0 e ) 

(11) 
(12) 

! sina, (13) 
(14) 



where 



eB 
2irm e c 

kT 



and K2 is the modified Bessel function of the second kind 
of the second order. Here, a is the angle between the 
wave vector of the emitted photon and the magnetic field. 
For our estimate, we use the average < sin a >= it/ 4. 

We assume that optical paths follow at least one loop 
along the circular photon orbit of a Schwarzschild black 
hole located at radius r = 3r g . The emitted intensity of 
radiation is then given by the expression 

4m=67rr ff j l/ . (15) 

The observed intensity is related to the emitted intensity 
by the third power of the redshift factor, which we take 
to be the gravitational redshift of a photon observed at 
infinity emitted from r — 3r g (i.e., we neglect for this 
simple estimate the high velocity of the flow, which will 
serve to increase the intensity at infinity). This yields: 



where 



lobs = (1 + z) 3 7cm, 

1 



l + z = 



= V / 3. 
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The photon ring has an approximate width of 
AL ~ 0.1r g = 0.51 /iarcsec 



(16) 
(17) 

(18) 
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Fig. 2. — Estimate of the flux density of the photon ring of 
Sgr A* assuming emission from a geometrically thin ADAF around 
a Schwarzschild black hole (curve labeled "model"). At wave- 
lengths A < 0.8 mm, the modeled flux density is of the order of 
0.2 Jy. We also plot (curve labeled "blackbody") a blackbody 
function at the same temperature and emission radius, which sets 
an upper flux density limit to the ring emission of the ADAF model 
at wavelengths A > 0.8 mm. Beyond this wavelength, the emission 
becomes optically thick. The solid line marks our estimate of the 
ring flux density in both wavelength ranges corresponding to the 
minimum of both functions. 

in the image plane (Johannsen & Psaltis 2010). Then, 
the subtended solid angle is given by the expression 



CI 



Lx AL 
D 2 



n 2 



(19) 



This yields our estimate for the observed flux density 
of the photon ring: 



^i/,ring — *ahs X CI. 



(20) 

In Figure we plot the modeled ring flux density as a 
function of wavelength. We also plot a blackbody func- 
tion evaluated at the same emission radius and temper- 
ature given by expression (J7J as an upper flux density 
limit. The above estimate exceeds the blackbody flux 
density at wavelengths A > 0.8 mm, and, thus, at these 
longer wavelengths, the emission becomes self-absorbed. 
Therefore, we use the minimum of these two flux densities 
as an estimate of the ring flux density. At wavelengths 
A < 0.8 mm, the modeled flux density of the photon ring 
of Sgr A* is ~ 0.2 Jy, about 1/15 of the total source flux 
density. Since A < 0.8 mm is also in the regime where 
VLBI observations are resolution-limited (see Figure [IJ , 
this range of wavelengths is optimal for such measure- 
ments. 

The ring diameter is determined solely by the mass 
via equation ([T]), and we can relate the black- hole mass 
and the angular diameter of the ring according to the 
expression 



Mo 



10 6 M m 



9.8 x 10" 



D 



/xarcsec kpc 



(21) 



Therefore, for VLBI imaging, the mass is proportional 
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Fig. 3. — 68% and 95% confidence contours of the mass and distance of Sgr A* for the combined distribution of stellar orbits and 
simulated VLBI measurements for a thermal noise-limited array at several wavelengths. Compared to (solid line) the 95% confidence 
contour of dynamical observations alone (Gillessen et al. 2009), the improvement of the mass and distance measurements is similar at all 
three wavelengths. For comparison, we also plot (left panel, dashed lines) the constant ratios M/D and M/D 2 in order to illustrate the 
dependence of both individual methods on the correleation between mass and distance. 



to the distance, and we can reduce the correlation be- 
tween mass and distance in combination with dynamical 
measurements . 

We now explore the prospect of combining dynamical 
measurements of Sgr A* with VLBI imaging observations 
of the photon ring. We analyze the best-case scenario of 
a thermal noise-limited VLBI array in order to assess 
whether such a measurement is worthwhile. Systematic 
limitations will degrade the VLBI measurement some- 
what, as discussed below. 

We employ a Bayesian analysis to estimate the proba- 
bility distribution over the mass and distance of Sgr A* 
from a measurement of the angular diameter of the pho- 
ton ring in combination with the constraints obtained 
from stellar dynamics. We take the latter as our prior, 
-Pprior(-Wj D), by converting into a likelihood the % 2 dis- 
tribution with v — 114 degrees of freedom obtained from 
the existing data set of the ephemerides of several S-stars 
(Gillessen et al. 2009). We assume a Gaussian posterior 
likelihood of obtaining a particular measurement of the 
angular diameter as 



PvLBi(data|M, D) 



1 



V2i 



■. exp 



%) 2 



2a 2 



(22) 



where # r ing is given by equation ([3]) and 9q = 53 /larcsec. 
We then use Bayes' theorem to write the likelihood of a 
particular mass and distance of the black hole given the 
data as 



P(M, £>|data) = CP V LBi(data|M,D)P prior (M, J D), 

(23) 

where C is the appropriate constant that normalizes the 
likelihood. 

The measurement uncertainty a is the key parameter 
of the likelihood Pvlbi- In the following, we estimate the 
scaling of the measurement uncertainty with resolution. 
The spatial frequencies of interest are those beyond the 
first null of the Bessel function describing the Fourier 
transform of the ring structure in visibility space. For a 
given resolution and ring size, there are N accessible half 



periods of the oscillation, where N is given by 

53 



N 



1 



21A E 



1. 



(24) 



The typical amplitude of a Bessel function oscillation of 
the first few maxima is of order 0.3 of the peak. There- 
fore, the signal we seek to measure has an amplitude of 



A = 0.3 x 



F v ,t 



nig 



xiV = 0.3x-xJV = 0.02iV, (25) 
3 



where inflow — 3 Jy is the observed flux density of Sgr A* 
near 1 mm (see, e.g., Broderick et al. 2009). 

We estimate the expected signal to noise ratio on this 
measurement by scaling it from the current VLBI mea- 
surements. Fish et al. (2011) measure the size of Sgr A* 
with a signal to noise ratio of ~ 40. Near-future VLBI 
arrays will incorporate 5 or 6 stations. One of them, 
ALMA, will have the sensitivity of 50 of the current sta- 
tions resulting in an overall increase in array sensitivity of 
a factor of up to roughly 9.3. In addition, the scheduled 
increase of recording bandwidth will increase the sensi- 
tivity by a factor of v8- To account for the variation in 
system temperature for typical observing conditions and 
receiver performance we introduce a degradation in SNR 
proportional to the observing wavelength squared, nor- 
malized to the 1.3 mm performance in Fish et al. (2011). 
The total signal to noise ratio of this measurement will 
then be 



SNR = 9.3 x V8 x 40 x A x 
/ 53 



A 



1.3 mm 



12 x A: 



V21A E 



1 



(26) 



with the width of the distribution given by the expression 
a - TTTr ~ 4.3 x A~l | 777-7— 1 ) /iarcsec. (27) 



SNR 



21A E 



In Figure [3l we plot confidence contours of the joint 
probability distribution for the combination of future 



MASSES OF NEARBY SUPERMASSIVE BLACK HOLES WITH VLBI 



5 



thermal noise-limited VLBI and current astrometric ob- 
servations at three different wavelengths. The solid line 
marks the 95% confidence contour determined by the 
stellar ephemerides (Gillessen et al. 2009). A VLBI 
measurement at a wavelength of 1 mm significantly im- 
proves the result from stellar orbits alone. At smaller 
wavelengths, the constraints on the mass and distance of 
Sgr A* are similar. In the rightmost panel of Figure 3, 
we have extrapolated the distribution width a given by 
expression (|27j) to a nominal wavelength of A = 0.5 mm. 
Measurements at such short wavelengths will be limited 
by weather conditions and may have to rely on a smaller 
array with fewer telescopes. 

Real observations will face more stringent limitations 
than those imposed by the interferometer thermal noise 
due to the complications of astrophysics and measure- 
ment systematics. The chief astrophysical limitation is 
the separation of the ring emission from the source struc- 
ture in the uv-plane. In our estimate, we have used the 
location of the nulls as a benchmark for the uncertain- 
ties we expect from the VLBI measurement. In practice, 
however, the full visibility function has to be analyzed 
with a pattern matching technique that identifies the 
structure of the ring. Such a technique has to extract 
the ring from a uv-plane that is only partially sampled 
by a given set of baselines. 

The physics of the accretion flow will also complicate 
things, as the structure of Sgr A* may vary over the 
course of an observation. However, because the ring 
structure is persistent and only weakly altered by rapid 
changes in the accretion flow we expect that tempo- 
ral averaging of the visibilities across multiple observing 
epochs will diminish the importance of such changes. 

The VLBI measurement itself must surmount system- 
atic limitations to make the moderate dynamic range 
measurements proposed here. Chief among these is the 
difficulty of calibrating the noise level at individual sta- 
tions, which imposed a 5% uncertainty in Fish et al. 
(2011) with the three-station array. In observations with 
the larger array considered here, there will be many more 
internal cross-checks available to improve the relative cal- 
ibration of stations (the absolute calibration is not im- 
portant). In particular, the use of three phased inter- 
ferometers (Mauna Kea, CARMA, ALMA) that simul- 
taneously record conventional interferometric data will 
permit scan-by-scan cross calibration of the amplitude 
scale of the array. 

Furthermore, the larger arrays will be able to make use 
of closure phases and closure amplitudes that are immune 
to calibration errors as part of the ring detection, al- 
though we have ignored such procedures here because of 
the difficulty of simply parameterizing the improvement 
they can permit. Other effects, such as the coherence of 
the reference systems between stations (reported as < 5% 
in Fish et al. 2011) can be more carefully measured and 
corrected to prevent them from imposing fundamental 
limitations to the ring detection. 

3. OTHER SOURCES 

Besides Sgr A*, there exist other nearby supermas- 
sive black holes, for which a combination of dynami- 
cal measurements and VLBI observations could be feasi- 
ble. Since these supermassive black holes are located in 
host galaxies other than the Milky Way, observations are 



much less affected by interstellar scattering. As an exam- 
ple, Broderick & Loeb (2009) and Takahashi & Mineshige 
(2011) analyzed the prospects of imaging the shadow of 
M87 with VLBI observations at several different wave- 
lengths. 
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Fig. 4. — Angular diameters and distances of several supermas- 
sive black holes. Sgr A* has the largest angular diameter, closely 
followed by M87 due to its high mass, making these sources ideal 
targets for VLBI imaging. Data taken from Giiltekin et al. (2009). 

In Figure 2J we plot the angular diameter of the pho- 
ton rings against the distances of a collection of nearby 
supermassive black holes. Sgr A* is closest to us and has 
the largest angular diameter, closely followed by M87 
and M31 due to their large black hole masses. The top 
dashed line indicates the resolution of a telescope array 
with a baseline equal to the diameter of the Earth (from 
equation ([5])) at a wavelength of 1 mm. For comparison, 
we also show the resolution of a future space telescope 
located at 30% the distance to the moon (comparable to 
the orbit of Chandra) at a wavelength of 0.5 mm corre- 
sponding to an angular diameter of about 1 /^arcsec as 
well as of the proposed Millimetron mission (Wild et al. 
2009) at a distance of 3 x 10 5 km and a wavelength of 
A = 0.4 mm. 

In order to be able to resolve the photon ring with 
VLBI, the key question is whether the accretion flow of 
the target supermassive black hole is optically thin. In 
some cases, the spectra of these sources peak at wave- 
lengths near A max ~ 1 mm, similar to the spectrum 
of Sgr A*, suggesting that the emission comes from an 
ADAF (Di Matteo et al. 2000; Doi et al. 2005). Naively, 
we would expect an approximately linear scaling of the 
electron density of an ADAF with the ratio M/M of the 
black hole, where M is its mass accretion rate. The de- 
tails of such a relation depend on a variety of factors, such 
as the temperature profile, the emissivity, and the radia- 
tive efficiency. However, most of the nearby supermassive 
black holes have very low radiative efficiencies (Ho 2009). 
Therefore, it is plausible that the accretion flows of these 
nearby supermassive black holes become optically thin at 
wavelengths that are comparable to 1 mm making them 
accessible to VLBI observations. Such observations are 
best carried out at wavelengths near the flux peak, where 
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Sources for VLBI Observations 
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et al. (2005); For S„: (1) VLBI: 1.3 mm (Doeleman et al. 2008); (2) VLBI: 86 GHz (Leo ct al. 2008); (3) Power law fit of mm-data from 
NASA Extragalactic Database, evaluated at 1 mm; (4) Submillimetrc Common-User Bolometer Array (SCUBA): 2 mm (Di Matteo et al. 1999); 
(5) Australia Telescope Compact Array (ATCA): 95 GHz (Sadler et al. 2008); (6) NObeyama Bolometer Array (NOBA): 2 mm (Leeuw et al. 2004); 
(7) Swedish-ESO 15m Submillimeter Telescope (SEST): 2 mm (Israel et al. 2008); (8) SCUBA: 0.85 mm (Bendo et al. 2006); (9) Very Long Baseline 
Array (VLBA): 8.4 GHz (Pellegrini ct al. 2003); (10) Plateau de Bure Interferometer (PdBI): 241.4 GHz (Schodcl ct al. 2007); (11) VLBI: 86 GHz 
(Middclberg ct al. 2005). 



the accretion flow is becoming optically thin. At wave- 
lengths A <C A max , the emitted flux is likely to be too low 
to be detected with a VLBI array, while at wavelengths 
A ^> A max , the accretion flow is optically thick. 

In Table 1, we summarize the angular diameters, dis- 
tances, masses, radio luminosities Lr, flux densities S u 
near 1 mm, ratios of the bolometric luminosity to the 
Eddington luminosity iboi/^Edd, and peak wavelengths 
Amax for supermassive black holes, whose photon rings 
have an angular diameter of at least 1 /zarcsec. 

In addition to Sgr A*, the black holes in the centers of 
M87, M31, and M60 are good potential targets for VLBI 
observations, because of the large angular diameters of 
their respective photon rings and, in the case of M87 
and M60, the measured peak in the synchrotron part of 
their spectra near 1 mm. M87 has a high measured flux 
density at 86 GHz (Lee et al. 2008) and should be readily 
observable at wavelengths close to 1 mm. In the case of 
M60, however, Di Matteo et al. (1999) report an upper 
limit on the flux density of 4 mJy. No similar flux density 
measurement of M31 has been reported to date. We 
estimate the flux density at 1 mm of M31 from a simple 
power law fit of mm-data from the NASA Extragalactic 
Database. This flux density is relatively small, and M31 
as well as M60 may be too faint to be observable with 
VLBI. Other sources, such as Centaurus A, are luminous 
enough to be detectable at wavelengths near 1 mm. With 
increasing VLBI resolution, even their photon rings may 
become observable. 

In the following, we assess the improvement on the 
mass measurements of the two supermassive black holes 
(Sgr A* and M87) whose photon rings have the largest 



angular diameters. We assume fixed distances of Dq = 
8.3 kpc (Gillessen et al. 2009) and D = 17 Mpc (Gultekin 
et al. 2009) for Sgr A* and M87, respectively. For 
Sgr A*, we estimate an error of the combined mass mea- 
surement from the existing stellar ephemerides and our 
simulated VLBI data of only ~ 5% (see Figure [3]). 

For M87, we estimate the smallest relative error that 
thermal noise-limited VLBI imaging observations of the 
ring can achieve from the signal to noise ratio for obser- 
vations of Sgr A* given by equation (f2"6"| . which we scale 
with the angular diameter of the black hole to obtain the 
expression 



SM 



/ rmg 



1 x-2 ( 



'ring 



12 



-A; 



V.21A n 



1 



(28) 



Note, however, that M87 has a much larger mass and that 
the dynamical timescales of its accretion flow is much 
longer. Therefore, VLBI imaging observations of its pho- 
ton ring will be much less affected by the variability of 
the accretion flow as in the case of Sgr A*. 

In Table 2, we compare the relative errors of the mass 
measurements of Sgr A* and M87 (Gillessen et al. 2009; 
Gultekin et al. 2009) with our estimate of the error of 
VLBI observations of the respective photon rings at sev- 
eral different wavelengths. In the case of Sgr A*, imaging 
its photon ring improves the error by a factor of about 
two. In the case of M87, imaging the photon ring at 
a wavelength of 0.5 mm would lead to a result that is 
similar to the current mass measurement. 

As we pointed out in section 2 for the case of 
Sgr A*, these errors require further refinement by in- 
depth imaging simulations. In addition, the morphology 
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Table 2. Relative Errors of Mass Measurements 



Source 


SM 1 


SM 1 


SM 1 


SM 


M | dyn 


M IvLBI 


M IvLBI 


M IVLBI 






(1.0 mm) 


(0.8 mm) 


(0.5 mm) 


Sgr A* 


0.09 


0.05 


0.05 


0.05 


M87 


0.29 


1.75 


0.42 


0.30 



of (sub-)mm VLBI emission can be complicated by the 
presence of a jet (for M87, see Broderick & Loeb 2009; 
Dexter et al. 2011). 

As in the case of Sgr A*, the combination of the results 
from both the dynamical and VLBI imaging observations 
of M87 would further reduce the error in the masses. The 
relative errors of mass measurements of both techniques 
likewise depend on the error in the measured distances 
to these sources. These errors, in turn, depend on un- 
certainties in the Hubble constant, peculiar motions of 
the gas in host galaxies, as well as assumptions on the 
proper motion of the Milky Way (see, e.g., Hodge 1981; 
Jacoby et al. 1992). The details of these effects on the 
relative errors of the mass are beyond the scope of our 
analysis. 

4. DISCUSSION 

In this paper, we investigated the prospects of mea- 
suring the mass and distance of Sgr A* as well as the 
masses of several other nearby supermassive black holes 
with a combination of dynamical observations and VLBI 
imaging of the respective photon rings of these sources. 
In order to resolve the photon ring of a black hole, its ac- 
cretion flow must be optically thin. We argued that the 
wavelengths at which the accretion flows of these sources 
become optically thin should be roughly comparable to 
the location of the peak in the synchrotron emission of 
Sgr A* and identified several supermassive black holes as 
optimal targets. 

We explored the prospects of imaging the photon ring 
of Sgr A* as well as of other nearby supermassive black 
holes with near-future VLBI arrays. We estimated the 
signal to noise ratio with which such arrays can image 
the photon ring in the best-case scenario if the VLBI ob- 
servations are limited by thermal noise. Based on our 
estimate, we simulated confidence contours of a mass 
measurement of Sgr A* using existing data of stellar 
ephemerides. We showed that the combination of both 
techniques can indeed reduce the correlation between 
mass and distance significantly resulting in relative er- 
rors of the mass and distance of only a few percent. We 
also identified several sources of uncertainty that have 
to be taken into account for an actual detection of the 
photon ring of Sgr A*. 



The uncertainties of measurements based on stellar or- 
bits will be further reduced by the continued monitoring 
and by the expected improvement in astrometry possi- 
ble with the second generation instrument GRAVITY Lor 
the Very Large Telescope Interferometer (Eisenhauer et 
al. 2011). Further improvements of the VLBI sensitiv- 
ity will be achieved by the Event Horizon Telescope, a 
planned global array of (sub-)mm telescopes (Doeleman 
et al. 2009a, b; Fish et al. 2009). 

We estimated the improvement of the mass measure- 
ment of M87 using VLBI techniques. Such observations 
are promising at wavelengths near 0.5 mm because of 
the large size of its photon ring. For M31 and M60, the 
supermassive black holes with the largest photon rings 
besides Sgr A* and M87, the flux densities may be too 
low to be detectable with VLBI. As the resolution of 
VLBI arrays increases, additional sources will become 
observable. 

Angular resolution of ~ 1 /larcsec requires longer base- 
lines and/or shorter observing wavelengths. However, 
the atmosphere precludes regular VLBI observations at 
wavelengths shorter than ~ 0.3 mm at even the best 
sites. 

A measurement of the photon-ring diameter is likely to 
yield useful results only if the observations extend well 
beyond the first null in the ring's visibility function. Even 
at a wavelength of A ~ 0.4 mm, space- VLBI observations 
will be required to reach this point for all except the first 
four entries in Table 1. 

VLBI between Earthbound antennas and a satellite has 
been achieved at 6 cm wavelength using the Japanese 
HALCA satellite (Hirabayashi et al. 1998). The recently 
launched RadioAstron (Kardashev 2009) will extend such 
observations to 1.2 cm wavelength and baselines as large 
as 4 x 10 5 km, for a resolution of < 10 /xarcsec. 

A future Explorer-class space mission designed to ob- 
serve at 1 mm wavelength or shorter, where source opac- 
ity and scattering effects will be far less important, could 
provide the angular resolution needed to study a far 
larger sample of sources. Such a capability may be 
provided by the Russian-European Millimetron mission, 
which plans to deploy a 12 m antenna with VLBI capa- 
bilities to 0.4 mm and maximum baseline > 3 x 10 5 km 
(Wild et al. 2009). 

We thank D. Zaritzky for useful comments. TJ and 
DP were supported by the NSF CAREER award NSF 
0746549. This work was supported at MIT Haystack Ob- 
servatory by NSF grants AST-0908731, AST-0905844, 
and AST-0807843. 
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